The initial conditions of the universe: how much isocurvature is allowed? 
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We investigate the constraints imposed by the current data on correlated mixtures of adiabatic and 
non-adiabatic primordial perturbations. We discover subtle flat directions in parameter space that 
tolerate large (~ 60%) contributions of non-adiabatic fluctuations. In particular, larger values of 
the baryon density and a spectral tilt are allowed. The cancellations in the degenerate directions 
are explored and the role of priors elucidated. 
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A major observational program is underway to deter- 
mine the physical properties of the universe. A key el- 
ement of this program is to fully characterise the pri- 
mordial fluctuations that seeded structure. The avail- 
ability of precision cosmic microwave background (CMB) 
datasets from satellite experiments [1] and smaller scale 
ground-based experiments [2], complemented by large- 
scale structure (LSS) data from galaxy redshift surveys 
[3,4], has now made it possible to work towards this am- 
bitious goal. 

The simplest characterisation of the primordial pertur- 
bations involves a Gaussian, nearly scale-invariant spec- 
trum of fluctuations arising from adiabatic initial con- 
ditions. These features are motivated by the simplest 
single-field models of inflation. To rigorously establish 
the character of the primordial perturbations, however, 
it is necessary to study a wider class of models. It there- 
fore seems reasonable to consider more general possibili- 
ties for the primordial perturbations. 

In this paper we adopt a phenomenological approach 
to the problem of determining the initial conditions of 
structure formation. Adiabatic initial conditions, char- 
acterised by a single, spatially uniform equation of state 
for the stress-energy content of the universe, have so far 
received much attention in the literature. Such fluctua- 
tions arise from perturbing all particle species in spatially 
uniform ratios, resulting in an overall curvature pertur- 
bation to the hypersurfaces of constant cosmic temper- 
ature. However, isocurvature perturbations also provide 
regular initial conditions for the evolution of fluctuations. 
Perturbations of this type result from spatially varying 
abundances of particle species that are arranged to cancel 
locally so that the curvature of the spatial hypersurface 
is unperturbed. 

In a universe filled with photons, neutrinos, baryons 
and a cold dark matter (CDM) component, four regu- 
lar isocurvature modes arise in addition to the familiar 
adiabatic mode. The cold dark matter isocurvature (CI) 
and the baryon isocurvature (BI) modes, in which spatial 
variations in the cold dark matter density or baryon den- 
sity are compensated by photon density perturbations, 
were studied some time ago [5] . The possibility of creat- 



ing a neutrino isocurvature density (NID) mode, in which 
perturbations in the neutrino density are balanced by op- 
posing photon density fluctuations, or a neutrino isocur- 
vature velocity (NIV) mode, where perturbations in the 
velocity of neutrinos are compensated by equal and op- 
posite photon velocity perturbations, was realised more 
recently [6,7]. 

When two or more perturbation modes are excited, 
the possibility of non-trivial correlations between these 
modes arises. The most general Gaussian perturbation 
of the five regular modes is completely characterised by 
the symmetric, matrix-valued power spectrum, 

P«(fc)-* 3 (k-k') = (i4<(k) A-(k')), 

where the indices (i, j — 1,2,3, 4, 5) label the modes and 
Ai are the mode amplitudes [6] . If attention is restricted 
to quadratic observables, like the CMB angular power- 
spectrum, this description also suffices for non-Gaussian 
fluctuations. Under the assumption that the power spec- 
trum for each mode is a smoothly varying function of k, 
each auto-correlation and cross-correlation mode may be 
parameterised by an amplitude and a spectral index. 

The viability of isocurvature perturbations in the light 
of observational data has been explored before [8,9], over 
restricted sets of parameters. A detailed exploration of 
the likelihood space of cosmological parameters and ini- 
tial conditions was undertaken in [10] to forecast how the 
WMAP and Planck satellite experiments would measure 
the amplitude of primordial isocurvature fluctuations. 
However, due to the unavailability of satellite data, this 
analysis was restricted to a perturbative exploration of 
the likelihood surface around an adiabatic model. Here 
we utilize current CMB anisotropy and large-scale struc- 
ture data to measure the amplitude of correlated non- 
adiabatic fluctuations over a broad range of flat cosmolo- 
gies. 

We consider a sixteen dimensional parameter space 
that includes the physical baryon and cold dark mat- 
ter density parameters, u>b and u> c , the fractional energy 
density of a cosmological constant, f^A, the optical depth 
parameter to the reionisation epoch, t, and 0, quantify- 
ing redshift-space distortions in the large-scale structure 
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FIG. 1. Marginalized distributions for the 'cosmological' 
quantities, for CMB (dashed), CMB+LSS (solid). Results for 
pure adiabatic models with CMB+LSS are dot-dashed. 



measurements. As a prior, we require r < 0.3, as in [11]. 
The mode correlations are described by the ten parame- 
ters, (Ai Aj), and a single overall tilt parameter, n s , for 
all the modes. To avoid introducing an artificial degen- 
eracy, we exclude the BI mode, as its prediction for the 
CMB spectrum is identical to that of the CI mode [10]. 

We treat the mode contributions as follows. Firstly, 
the CMB and matter power spectra for each pure mode 
are normalized by its mean square power, defined as its 
contribution, (ST 2 ) = *£(2t + l)C e , to the CMB tem- 
perature anisotropy from £ m i n — 2 to £ max = 2000. 
The cross-correlation spectra are rescaled by the geo- 
metric mean of the corresponding pure mode powers. A 
mixed model is constructed from these spectra using co- 
efficients, Zij, that are normalized to have total rms am- 
plitude equal to one. Geometrically a set of symmetric 
corresponds to a position z, on the unit nine-sphere. 
The CMB and matter power spectra for the mixed model 
are then rescaled by the mean square power of the admix- 
ture before being multiplied by a normalization parame- 
ter, P, to yield the final model spectra. The parameter P 
quantifies the mean square power of the resulting model, 
which is well measured from current data. We impose 
uniform priors on the direction z, and on P. 

The CMB temperature and temperature-polarisation 
cross-correlation angular power spectra, and the matter 
power spectrum were evaluated using a modified version 
of the grid-based DASh software package [12], extended 
to compute auto-correlation and cross-correlation spectra 
for adiabatic and isocurvature modes. The computation 
of 10 mode spectra, for a given cosmological model, takes 
20 seconds on a single Pentium 3 processor. The CMB 
power spectra computed with our improved version of 
DASh agree with those computed using the latest version 
of CMBFAST [13] to within 0.5%, which is sufficient for 
our likelihood computations. 

To sample the parameter space efficiently, we employ a 
Markov Chain Monte Carlo method, using the Metropolis 
algorithm. We are able to sample the posterior distribu- 
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FIG. 2. Top four panels: marginalized distributions of 
relative power contributions zu of the four auto-correlated 
modes: adiabatic (AD), cold dark matter density isocurva- 
ture (CI), neutrino velocity isocurvature (NIV) and neutrino 
density isocurvature (NID) for CMB (dashed) and CMB+LSS 
(solid) . Bottom six panels: marginalized distributions for rel- 
ative power contributions of cross-correlated modes, for 
CMB (dashed) and CMB+LSS (solid). 

tion for each parameter quickly and accurately using an 
efficient proposal distribution, for which we take a scaled 
best guess of the covariance matrix of the posterior. To 
obtain an adequate estimate of this covariance matrix, 
requires updating a series of chains, of cumulative length 
i=a 10 5 steps. However, with this approximate covariance 
matrix, the subsequent single chains used in our analysis 
converge over much shorter lengths, requiring typically 
20,000 likelihood calculations. Our criterion for conver- 
gence to the underlying distribution is that the sample 
variance of the mean of each parameter, computed us- 
ing a spectral method, is less than one percent of the 
variance of that parameter. More details of the above 
methods are given in [14]. 

We compare our model spectra to measurements of 
the CMB anisotropy from WMAP [1] and a compila- 
tion of CMB data on small scales [15] from the ACBAR, 
BOOMERANG, CBI, and VSA experiments [2]. A com- 
bination of the galaxy power spectrum measurement by 
2dFGRS [3] and CMB measurements constitutes our sec- 
ond dataset. We use the CMB likelihood functions pro- 
vided in [16] and [15] for the WMAP and small-scale 
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FIG. 3. The CMB temperature spectrum (panel 1) 
and the temperature-polarization cross-correlation spec- 
trum (panel 2) are plotted for the maximum likeli- 
hood adiabatic model (dashed line) with parameter values 
(cut, La c ,n A ,n a ,T,(3)=(0.023, 0.12, 0.72, 0.97,0.14, 0.49) and a 
high-likelihood mixed model (solid line) with parameter val- 
ues (0.041,0.13,0.75,1.06,0.28,0.37), with the CMB data 
overplotted. The galaxy power spectrum P(k) for these mod- 
els along with the pure mode contributions to the mixed 
model is shown in panel 3, including sizes of LSS error bars. 

CMB data, respectively, and the large-scale structure 
likelihood function available in [3] , restricting ourselves to 
the linear regime 0.01 < k < 0.15 /iMpc -1 . The normal- 
ization of the observed galaxy power spectrum depends 
on the rcdshift-space distortion parameter, /?, that was 
measured in [17] to be 0.43 ± 0.07, which we include as 
a Gaussian prior. We follow the method of [16] in nor- 
malizing the theoretical matter power spectrum to the 
galaxy power spectrum. However, we do not include an 
independent prior on the bias, 6, nor treat the effects of 
peculiar velocities, important at smaller scales. 

For pure adiabatic models our results are consistent 
with the WMAP analysis [11], for CMB data alone and 
for CMB data combined with LSS data. The quality of 
fit of the highest likelihood models, together with the 
median parameter values and 68% confidence intervals 
are given in table I. The marginalised parameter distri- 
butions are also displayed in Fig. 1. When correlated 
adiabatic and isocurvature initial conditions are allowed 
the situation changes significantly. We find that large 
fractions of isocurvature models are tolerated by current 
data. The best-fit mixed model has a likelihood higher 
by a factor f* e 5 / 2 than that of the best-fit pure adiabatic 
model. Fig. 3 shows the best-fit model spectra with their 
parameter values listed in the caption; the pure adiabatic 
model and the mixed model are essentially indistinguish- 
able. The relative magnitude of the non-ad iabatic com- 
ponent z lso , defined as z lso = — zJ D ), is measured 
for the CMB data alone to be 0.84+q"| and comprises 
roughly equal proportions of pure isocurvature modes. 
To extract a fractional non-adiabatic contribution we de- 
fine / ISO = z ISO /(z ISO + z A d), calculated to be 0.60t ; n . 

If large-scale structure data is included the relative 
magnitude of non-adiabatic modes is slightly reduced to 
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TABLE I. Median parameter values, 68% confidence inter- 
vals, and best fit \ 2 P er degree of freedom v, for pure adiabatic 
models (columns 1 and 2) and mixed models (columns 3 and 
4). The coefficients Zij are tabulated for the different modes. 

the median value of 0.79±o;°3 , with / ISO = 0.57 ± 0.09. 
The inclusion of isocurvature modes does not improve the 
fit to the 2dF data, but neither does this dataset signif- 
icantly constrain the isocurvature fraction. This is true 
because we have allowed the normalization of the galaxy 
power spectrum to vary, only indirectly constraining it 
through our prior on j3. We find that smaller values of (3 
are favoured, corresponding to larger biases at fixed Q m . 
This is not surprising given that the matter power spectra 
generated in isocurvature models are lower in amplitude 
relative to the adiabatic spectrum, as illustrated in Fig. 
3. Therefore models with higher isocurvature content 
have correspondingly larger biases. For both datasets 
the median value of n s is larger than in the adiabatic 
case with a much broader distribution. 

Turning to the cosmological parameters, we find that 
lu c and f^A have broader distributions shifted relative 
to the adiabatic case, when CMB data alone is con- 
sidered. The constraints on these parameters become 
tighter when 2dF data is included, mainly because the 
matter power spectrum is sensitive to a combination of 
fl m and h, the dimcnsionless Hubble parameter, which 
in our parameterisation translates into improved con- 
straints on ijj c and within the range of flat mod- 
els. The reionisation parameter, r, which has a fairly 
broad distribution in the adiabatic case, is very poorly 
constrained in mixed models. 

The greatest impact, however, is on the baryon den- 
sity distribution, which is significantly broadened, with 
a median value twice as large as in the adiabatic case. 
The drastic shift in the baryon density is associated with 
a flat direction in parameter space, which was observed 
previously in [9] but not investigated further. Here we 
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FIG. 4. Top panel: The total derivative along the degener- 
ate direction (dot-dashed) described in the text is decomposed 
into variations in four parameters (a;;,, n a , AD, NIV), which 
cancel to within the experimental error bars (shaded). Bot- 
tom panel: (left) marginalized prior for the adiabatic mode 
contribution z A d, no data; (centre) marginalized posterior for 
zad with uniform prior on all models (dashed) and with uni- 
form prior on z A d (solid); (right) marginalized posterior for 
cot with priors as in centre panel. 

clarify the nature of this degeneracy. In Fig. 4 we illus- 
trate how a change in a high-likelihood model Cg due to 
a 10% increase in Ub is compensated by a 40% increase in 
the power of the pure NIV mode, a 1.4% increase in the 
value of n s and a 14% decrease in the power of the pure 
AD mode, to well within the CMB error bars. We have 
found that there is a high-likelihood plateau connecting 
the mixed model to the adiabatic model. It is along this 
plateau that Ub and n s are shifted to higher values. 

We found that the allowed fractional contribution of 
isocurvature modes depends sensitively on the choice of 
prior. The problem is that no natural notion of a "uni- 
form" prior stands out, in particular when correlations 
between the modes are allowed. The prior presented ear- 
lier can be argued to be biased against models where 
a single mode dominates (such as the pure adiabatic 
model) because such models have a comparatively minus- 
cule phase space density. This "entropy" effect is illus- 
trated in the bottom left hand panel of Fig. 4, where the 
a posteriori distribution for z AD that would result with 
no data is plotted. We observe that the pure adiabatic 
model has zero density solely on account of phase space 
volume effects. In the bottom centre panel of Fig. 4 we 
show the posterior distribution for z AD that would result 
from dividing out this phase space density effect (solid 
curve), compared to the original posterior (dashed). Af- 
ter the rescaling of the prior, the allowed isocurvature 
contribution is greatly suppressed. In the bottom right 
panel of Fig. 4 we see that the reweighted posterior for 
u>b overlaps with the conventional value, although much 
higher values are still allowed. 

We have restricted our analysis to flat models and have 
assigned a single spectral index to all modes. Relaxing 
either of these restrictions will open up interesting re- 
gions of parameter space. However, it is evident from 



the results presented here that further precision data 
is required before meaningful constraints can be placed 
on extended parameter sets. Future high quality CMB 
polarisation data promises to provide such strong con- 
straints [18], while improved temperature and polarisa- 
tion spectra from WMAP, the recently available SDSS 
dataset and higher quality small-scale temperature mea- 
surements from ground-based experiments will allow fur- 
ther progress on these issues. 
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